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ABSTRACT 

The optical-IR images of the Northern and Southern Hubble Deep Fields are used to measure 
the spectral and density evolution of early-type galaxies. The mean optical SED is found to 
evolve passively towards a mid F-star dominated spectrum by z~2. We demonstrate with real- 
istic simulations that hotter ellipticals would be readily visible if evolution progressed blueward 
and brightward at z>2, following a standard IMF. The colour distributions are best fitted by 
a 'red' IMF, deficient above ~ 2Mq and with a spread of formation in the range 1.5< Zf <2.5. 
Traditional age dating is spurious in this context, a distant elliptical can be young but appear 
red, with an apparent age > 3Gyrs independent of its formation redshift. Regarding density 
evolution, we demonstrate that the sharp decline in numbers claimed at z>l results from a 
selection bias against distant red galaxies in the optical, where the fiux is too weak for morpho- 
logical classification, but is remedied with relatively modest IR exposures revealing a roughly 
constant space density to z~2. We point out that the lack of high mass star-formation inferred 
here and the requirement of metals implicates cooling-flows of pre-enriched gas in the creation 
of the stellar content of spheroidal galaxies. Deep-fleld X-ray images will be very helpful to 
examine this possibility. 

Subject headings: cosmology: elliptical galaxies — cosmology: observations — galaxies: 
starformation and redshifts 



port for early dust is controversial. Claims of high 
star-formation rates from Far-IR imaging of the HDF 
(Hughes et al 1998) have been corrected in the ra- 
dio (Richards 1998) and identified with low redshift 
disk galaxies and AGN. This is consistent with the 
non-detection of Far-IR emission from targeted ob- 
servations of optically luminous z ~ 3 galaxies (Scott 
etal. 1998). 

In the context of hierarchical models, it is natu- 
ral to view E/SO galaxies as the end product of a 
merging process and hence to predict declining num- 
bers with increasing redshift. Locally at least, merg- 
ing of disk galaxies is seen to create some spheroidal 
shaped objects (Schweizer & Seitzer 1992). At faint 
magnitudes claims have been made for a decline in 
the space density of red and/or elliptical galaxies at 
z>l, (Kauffmann, White, & Chariot 1996; Zepf 1997; 
Franceschini et al. 1997; Kauffmann & Chariot 1998; 
Menanteau etal. 1999; Barger et al. 1999) mainly 

^Based on observations made with NASA/ESA Hubble Space Telescope which is operated by AURA, Inc., under contract with 
NASA. 

^ Based on observations made at the Kitt Peak National Observatory, National Optical Astronomy Observatories, which is 
operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation. 

''Observations have been carried out using the ESO New Technology Telescope (NTT) at the La Silla observatory under Program- 
ID Nos. 59.A-9005(A) and 60.A-9005(A). 



1. INTRODUCTION 

Elliptical galaxies are anomalous in many respects 
when considered in the context of the standard 
ideas regarding galaxy and star-formation. De- 
spite the absence of star-formation today only min- 
imal passive evolution has been identified to z~l, 
mainly from optical-IR colours of cluster sequences, 
which are marginally bluer than k-corrections pre- 
dict (Stanford, Eisenhardt & Dickinson 1998). At 
higher redshift, examples of luminous red galaxies are 
found with F and G-star dominated spectra (Dunlop 
etal. 1996; Spinrad etal. 1998; Broadhurst & Frye 
1999). No bright blue precursors have been identi- 
fied. The absence of precursors naively implies the 
early epoch was obscured by dust or restricted to 
unobservably high redshifts. Benitez etal. (1999) 
strongly limit any unobscured formation to z>10, in 
the deepest available VLT/NICMOS images. Sup- 
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on the basis of optical imaging. Conflicting with this 
is the simplicity of structural and color relations be- 
tween elliptical galaxies (Fabcr &: Jackson 1976), par- 
ticularly in rich clusters where coeval monolithic for- 
mation is inferred (Bower, Lucey, & Ellis 1992). 

If gas rich mergers of spirals are to produce the 
ellipticals then the enhanced alpha-element abun- 
dance generated by a brief merger induced episode 
of star-formation is unacceptably diluted by the pre- 
existing SNIa-enriched ISM (Thomas, Greggio & 
Bender 1999). Furthermore the mass of stars formed 
during a merger is limited by the general absence of 
an intermediate-age stellar population in post-merger 
ellipticals (Silva &: Bothun 1998; James & Mobasher 
1999). Metallicity is also a problem with monolithic 
collapse formation, as closed box star-formation does 
not account for the observed lack of low metallicity 
stars (Worthey 1994) but implies pre-enrichment of 
the gas (Thomas, Greggio & Bender 1999). 

Recently, distant red ellipticals and other spheroid 
dominated galaxies at l<z<2 have been detected in 
the deepest combination of optical-IR images - a 
smaU NICMOS/VLT field (Benitez etal. 1999, TVeu 
etal. 1998, Stiavelli etal. 1998). Here we analyse a 
much larger sample of distant red galaxies by com- 
bining optical-IR photometry from both Hubble Deep 
Fields, to measure the rates of spectral and density 
evolution (§3, §4) with photometric redshift measure- 
ments (§2) and discuss new implications for their for- 
mation (§4). 

2. OBSERVATIONS 

The observations used here arc the deep HST op- 
tical images in the UBVI (Williams et al. 1996, 
1999) and the JHK images from KPNO (Dickinson 
etal. 1997) and from the SOFI instrument on the 
NTT in the south (da Costa etal. 1998). We use 
the published zero-points, filter transmission and de- 
tector response curves. Magnitudes are measured us- 
ing SExtractor (Bertin & Arnouts 1996), and pho- 
tometric redshifts are estimated by maximizing the 
likelihood with respect to a set of instantaneous 
burst spectra calculated using the [Z/Zq] = —0.2 
Bruzual/ Chariot spectrophotometric package (Lei- 
therer et al. 1996). Stars are distinguishable to very 
faint magnitudes by both a stcllarity index (Bertin & 
Arnouts 1996) and, interestingly for red stars, by a 
poor fit to redshifted red galaxy spectra. 

3. SPECTRAL EVOLUTION 

The first point to note is that the choice of metallic- 
ity does not significantly affect the redshift estimate. 
The break at 4000A is so sharp that with accurate 
optical-IR magnitudes the redshift can be determined 
to ~15% with a ruler. A comparison of photo-z's with 



the 10 spectroscopic redshifts of elliptical galaxies the 
HDF is very good with a hint of a Sz=0.1 systematic 
overestimate. 48 objects are found with SEDs well 
matched to early type galaxies. The majority clearly 
fit a de-Vaucouleurs profiles and not an exponential 
(Fig 2). The remainder are too faint in the optical 
and of too low resolution in the IR to constrain the 
morphology. In other words, our sample contains all 
objects consistent with a de-Vaucoulours profile and 
a passively evolving SED. Half of these galaxies lie 
at z>l and extend to z=2.5. Few if any moderately 
bright blue ellipticals (~ 5 at mj ~ 23 — 24) are 
missed this way. 

The spectra are compared in the restframe (Fig 1). 
A clear evolutionary trend emerges towards a mid F- 
star dominated spectrum by z '~ 2. Hotter A-star 
dominated spectra would be very easily recognized 
if ellipticals evolved further according to a standard 
IMF at z>2, such young (<lGyr) precursors being 
very bright and blue. This simple result suggests 
that the passive evolution of elliptical galaxies be- 
gins at ~ 2Mq. For approximately ~l-2Gyrs after 
formation the spectrum of such a stellar population 
has no detectable spectral evolution, consistent with 
the slow evolution found here above z ^ 1 (Fig 1). 
This level of evolution corresponds to a change in 
magnitude in the rest-frame B of only l.™2 between 
z~2 and the present. A small but detectable variance 
among SEDs is observed at any redshift (Fig 1) with 
evidence for greater homogeneity at high redshift. 

The bluest 3 objects marking the starting point of 
the color-colour tracks in Fig 3 arc estimated to lie 
in the range l<z<2. These objects contain a small 
blue excess in U and B relative to an F-star spectrum 
(Fig 2, bottom panel) which is spatially distributed 
like the general light profile (Fig 2), ruling out an 
AGN contribution. Accommodating this with some 
A-star light steepens the IR appreciably, requiring a 
redder IMF for a good fit. Nebular continuum emis- 
sion is an interesting possibility. Spectroscopy would 
be very helpful in understanding these relatively blue 
ellipticals. 

4. DENSITY EVOLUTION 

A proper assessment of density evolution requires 
simulated images to account for the very strong red- 
shift dependent k-correction. Simulations are made 
in all bands using the local luminosity function of 
early type galaxies (Pozzetti et al 1998) and matched 
in background noise, pixel scale and PSF of each pass- 
band, using a variation of the machinery described in 
Bouwens, Broadhurst & Silk (1998). Selection and 
photometry of both the observed and simulated im- 
ages is performed identically. Fig 3 shows a compar- 
ison with a model in which the density is fixed and 
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only the observed minimal spectral evolution takes 
place to Zf=2.5. It is clear that the numbers and 
luminosities of red galaxies has not changed much 
between z~2 and the present, in agreement with the 
claim of Benitez etal. (1999), but inconsistent with 
other estimates, in particular previous optical work. 
The unknown volume at high redshift translates into 
an uncertainty in the predicted numbers at z>l, so 
that both low 0, and flat A dominated models under- 
predict the data by ~30%. 

5. DISCUSSION AND CONCLUSIONS 

Our findings show that the passive evolution of el- 
lipticals evolves slowly to a mid-F star spectrum by 
z ^ 2. Bluer ellipticals are conspicuous by their ab- 
sence, and at face value, this simply suggests that the 
main sequence in elliptical galaxies does not extend 
above ~ 2Mq. It is also clear that most elliptical 
galaxies form at z>l given the lack of any significant 
decline in their space density with redshift, subject to 
a factor of 30% uncertainty from the unknown vol- 
ume. These conclusions are surprising given the high 
metal content of ellipticals and implies some gas pre- 
enrichment. This requirement is more palatable in 
light of recent evidence of outflows in higher redshift 
star forming galaxies, in particular lensed galaxies for 
which there is sufficient signal to detect blueshifted 
ISM absorption lines (Franx etal. 1997; Frye & 
Broadhurst; Frye etal. 1999). Such outflows will 
be preferentially enriched with alpha-elements from 
SNII activity. 

Independent of the observed outflows, Renzini 



(1997) has argued convincingly that SNII enrichment 
of the ICM is indicated by the predominantly alpha- 
element enriched gas. In the context of hierarchical 
evolution early enriched material will cool onto the 
later forming more massive halos. Locally, examples 
of near solar enriched cooling-flow X-ray gas is found 
in groups and clusters of galaxies centered on giant 
ellipticals of similar metallicity (Finoguenov & Pon- 
man 1999). We suggest that cooling may be respon- 
sible for the formation of the stellar content of ellip- 
tical galaxies more generally, naturally leading to a 
bottom heavy IMF consistent with our results. Invit- 
ing this simple picture is the remarkable correspon- 
dence between the most luminous X-ray cooling gas 
with impressively large cD galaxies (Fabian, Nulsen, 
&: Canizares 1991). Hence it is perhaps not surprising 
to find that such objects contain young stellar popu- 
lations (Mehlert et al. 1997) if this cooling gas forms 
visible stars. A clear test of the possible role of cool- 
ing flows in the formation of spheroidal galaxies will 
be provided soon by deep field X-ray imaging, like 
the planned deep AXAF field (Giacconi etal. 1999). 
Constraining the numbers of even higher redshift red 
galaxies requires deeper IR imaging to explore beyond 
z=2. The ISAAC camera on the VLT has the area 
and efficiency to achieve this, extending ~ 2 magni- 
tudes fainter than the relatively low resolution ~4m 
IR imaging used here. 

We thank Piero Rosati and Alvio Renzini for use- 
ful conversations. TJB acknowledges NASA grant 
AR07522.01-96A. 



REFERENCES 



Barger, A.J., et al. 1998, ApJ, 501, 522. 

Barger, A. J., Cowie, L.L., Trentham, N., Fulton, E. Hu, E.M., 
Songalia, ^ A Hall^ T) 1999 ApJ, in press, submitted to 
MNRAS, 

Bemtez, N., 



astro- ph/980929c . 

^rn arlhiirst, T..I.,"Rniiw ens. R.J., Silk, J., & Rosati, 



R 1999, ApJ, |astro-ph/9811129 . 
Berlin, E., & Arnouts. S. 1996, A&AS, 117, 393. 
Bower, R.G., Lucey, J.R., & Ellis, R.S. 1992, MNRAS, 254, 601 
Bouwens, R.J., Broadhurst, T.J., & Silk, J. 1998, ApJ, 506, 

557. 

da Costa, F. Nonino, R., Re 



igelink, R. etal. 1998 A&A 

submitted, astro-ph/9812105 

Dickinson, M., et al. 1997, ^ttp:/ /archive. stsci.edu/hdf/hdfirim. html 



Dunlop, J., Peacock, J., Spinrad, H., Dey A., Jimenez, R, 
Stern, D., & Windhorst, R.A. 1996, Nature, 381, 581. 

Fabian, A.C., Nulsen RE. J., & Canizares, C.R. 1991, 
A&A Rev., v2, 191. 

Faber, S. & Jackson, R. 1976, ApJ, 204 668. , 

Finog uenov. A . & Ponman, T.J. 1999, MNRAS, in press, |astro- 
ph/ ]990110C| . 

Franceschini, A., Silva, L., Fasano, G., Granato, L., Bressan 
A., Arnouts, S., & Danese, L. 1998, ApJ, 506, 600. 

Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950. 

Franx, M., lUingworth, G.D., Kelson, D., VanDokkum, P.G. 
Pieter G., Tran, K. 1997, ApJ, 496, L75. 

Fruchter et al. 1999, to appear in AJ. 

Frye, B. & Broadhurst, T. 1998, ApJ, 499, L115. 

Hughes, D.H., etal. . 1998, Nature,394,241 



James, P. & Mobasher, B. 1999, astro-ph/9901401 
Kauffmann, G. & Chariot, S. 1998, MNRAS, 297, L23. 
Kauffmann, G., Chariot, S., & White, S.D.M. 1996, MNRAS, 

283, L117. 
Leitherer et al. 1996, PASP, 108, 996. 

Mehlert, D., Bender R., Saglia, S., & Wegner, G. 1997, A&A, 
332, 33. 

Menanteau, F., Ellis, R.S., Abraham, R.. G.. Barger. A.J.. fc 
Cowie, L.L. 1999, MNRAS, submitted, |astro-ph/9811465 . 

Pozzetti, L., Bruzual, G. & Zamorani, G. 1996, MNRAS, 'ikl, 
953. 

Spinrad, H., Dey, A, Stern, D., Dunlop, J., Peacock, J., Jiminez, 
R., & Windhorst, R. 1998, MNRAS, 296, 1089. 

Stanford, S.A., & Eisenhardt, P.R., & Dickinson, M. 1998, ApJ, 
492, 461. 

Renzini, A. 1997, ApJ, 488, 35. 

Richards, E.A. 1998, ii i Before the Dar k Ages: When the 
Galaxies were Young, |astro-ph/9811120 . ^ 

StiaveUi, M., et al. 1999, A&A, submitted, |astro-ph/9812102 . 

Schweizer F. & Seitzer. P. 1992, AJ, 104, 1039. 

Silva, D. R. & Bothun, G. D. 1998, AJ, 116, 85. 

Thomas, D., Greggio, L., & Bender, R. 1999, MNRAS, 302, 
537. 

Treu, T., et al. 1998, A&A, 340, LIO. 
Zepf, S. 1997, Nature, 390, 377. 
Williams, R., et al. 1996, AJ, 112, 1335. 
Williams, R., et al. 1999, in preparation. 
Worthey, G., 1994, ApJS, 95, 107. 



4 



T r 




5000 1.5x104 2x104 5000 1.5x10* 2x10* 

Wavelength (Angstroms) lambda (rest-frame) 

Fig. 1. — The spectral evolution of the early type galaxies can be traced with redshift via their mean SEDs (a-e) in the range 
z<2.5. A detectable variance between SEDs is apparent at all redshifts but a general smooth trend is apparent towards the 
peaky spectrum of an F-star dominated spectrum. Panel (f) shows the trend of the mean model SED over the full redshift range 
normalized in the red. Nothing bluer than a mid type F-star spectrum is found. 
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Fig. 2. — Typical spheroidal galaxies as a function of redshift. We include the optical colour image, the IR (KPNO or NTT) 
image, the SED fits (model fiuxes are given by open triangles) to the observations (square points), and the observed 1-D profiles 
(crosses with errors) compared to the best-fit de-Vaucouleur profile (solid line) and exponential profile (dotted line). The bottom 
panel shows an example of an anomalous morphologically early-type galaxy with a U-B excess (see text). 



6 




Fig. 3. The dot plot on the right compares the rcdshift and absolute magnitudes (at 5000 A) of the observations (solid 
circles) with ellipticals recovered from our simulations (open circles, 2.5x in density) for our best-fit passive evolution model 
truncated at 3Af0 (using the Fioc & Rocca-Volmerange (1997) spectrophotometric tables) with formation redshifts distributed 
between z=1.5 and 2.5. The histogram shows the observed redshift distribution compared with the above best-fit model (dashed 
line) and a passive evolution model with Zf = 3 for two choices of geometry (O = 1, thick line; f2 = 0.3, thin line). Clearly, at 
most a factor of 30% decline in red galaxies is measured for the large volume models, but no trend to lower luminosity is found 
with increasing redshift. Colour-colour diagrams are also shown with evolutionary tracks indicating the sensitivity to formation 
redshift, Zf =1.5,1.75 & 2, progressing redder for later formation. 



